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QCD Sum Rules on the Lattice 

Chris R. AUton'' * and Stefano Capitani*^ 
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We study the work of Leinweber by applying the Continuum Model of QCD Sum Rules (QCDSR) to the 
analysis of (quenched) lattice correlation functions. We expand upon his work in several areas and find that, 
while the QCDSR Continuum Model very adequately fits lattice data, it does so only for non-physical values of its 
parameters. The non-relativistic model is found to predict essentially the same form for the correlation functions 
as the QCDSR Continuum Model but without the latter's restrictions. By fitting lattice data to a general form 
which includes the non-relativistic quark model as a special case, we confirm it as the model of choice. 
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1. INTRODUCTION 

Lattice Gauge Theory and QCD Sum Rules 
(QCDSR) are two areas of research which have 
been widely employed to deepen the understand- 
ing of systems governed by the strong interac- 
tions. In this paper we combine both methods 
following the pioneering work of Leinweber . 
Utilising known results in QCD Sum Rules we 
apply them to Lattice Gauge Theory in order to 
check that the two are compatible. 

We extend the work of [|l]j2| by including 
mesonic states, by applying our analysis to sev- 
eral P values, and by studying data from both the 
Wilson and SW-Clover actions. We also study the 
non-relativistic quark model predictions. 

The lattice data used is from the APE collabo- 
ration and is a mixture of (quenched) Wilson and 
SW-Clover with beta values ranging from /? — 6.0 
to /3 = 6.4. Zero-momentum 2-point hadronic 
correlation functions were studied for the nucleon, 
delta and the following mesonic operators: P5, Aq 
& Vi (in an obvious notation). A longer version 
of this work appears in 

2. QCDSR CONTINUUM MODEL 

The quark propagator in Euclidean space's 
Wilson OPE reads (in the coordinate gauge), 
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This is substituted into the (Wick contracted) 2- 
point function leading to the generic form: 

Cr'^W- E 7^-C„0„(™„(:gg:)), (2) 
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where C„ is a numerical coefficient, 0„ is some 
operator function, and no is a positive integer. 
It is now useful to express the timesliced corre- 
lation fimction, 6*2 (i), in the spectral representa- 
tion: G2{t) = p{s)e-'*ds. The OPE spectral 
density, p^^^{s), is calculated by means of a sim- 
ple inverse Laplace transform of G^^^, 
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where Re{t) > 0. In the sum over n, we have 
kept only the positive values of n since we are 
interested in the leading terms as t — *■ 0. 

The QCD Continuum Model is introduced by 
setting a threshold sq in the energy scale s, so that 
the excited states' contribution to G2{t) is given 
only by the energies above that scale. Performing 
the integral over s we obtain. 
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The full correlation function contains the ground 
state as well as the above continuum contribution 
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Qcont^ This is simply included as a (5 function in 
p{s), leading to the full correlation function 

gF'^HO = ^^-'"+^ p°^^(.)e--d.. (5) 

There are four parameters in this Ansatz for 
G^^^^{t): Z and M for the ground state, sq (the 
continuum threshold), and ^. The parameter ^ 
is introduced to allow for lattice "distortions" in 
the normalisation of the excited states ^ . In the 
continuum limit ^ should be equal to one. 

We have calculated the OPE expansions for all 
the channels studied. These are listed in 

3. FITS TO DATA 

Correlation functions were fitted to the follow- 
ing three functional forms, F{t): 
1. Conventional Single State Fit ("1-exp") 

1. e. using F{t) = ^er^'K 

2. QCD Continuum Model Fit ("Cont") 
These are described in the previous section. 

3. Conventional Two State Fit ("2-exp") i.e. 
using F{t) = affC"*" + ^e^^'^'K This func- 
tional form was chosen since it is traditionally 
used as the generalisation of "1-exp" to include 
the higher mass state(s). It has the same number 
of parameters as the "Cont" fit and the results of 
these two fits will be directly compared. 

The "1-exp" fits are used here to gain "stan- 
dard values" for Z and M and are fitted in 
the asymptotic region only. In the case of the 
QCD Continuum Model fits, we use time win- 
dows which begin very close to the origin, e.g. 
t = 2 — 28 for /3 = 6.0. So that a direct compari- 
son can be performed, the same time windows are 
used to perform the "2-exp" fits. 

The full results of the fits to the three fitting 
functions for all the lattice correlation data are 
presented in [|| . We show in Table |l|, the fits for 
the vector meson channel for the Wilson /5 = 6.1 
only. (The other fits follow a similar pattern.) 

The main features of the fits are now sum- 
marised. The Continuum Model reproduces the 
data better than the 2-exp case, i.e. (i) its 
values are lower (sometimes by an order of mag- 
nitude), and, (ii) its ground state parameters Z 
and M are closer to the "1-exp" control case. 



F{t) 


1-exp 


Cont 


2-exp 


Z(xlO-^) 


0.57(3) 


0.60(2) 


0.68(2) 


Ma 


0.466(2) 


0.468(2) 


0.474(2) 


So [GeV] 




1.9(1) 








-3.25(4) 




xVd.o.f. 


0.4(6)/ll 


9(5)/22 


110(20)/22 



Table 1 



Values for the fitting parameters for the vector 
meson (Wilson f3 — 6.1 case). 



We also note that the 2-exp fit values over- 
estimate Z and M in every channel studied. In 
1^ we give an explanation for this effect. 

Clearly it is not sufficient for the "Cont" fits 
to reproduce the ground state parameters Z and 
M, and to have a sensible x^i they must also give 
values for the continuum parameters, sq and ^ 
which are acceptable within the assumptions of 
the QCDSR Continuum Model. These are: (a) 
since sq corresponds to a physical threshold, they 
should be constant (in GeV) for each channel as 
f3 is varied; (b) sq should be large enough to be in 
a region where perturbation theory is valid; and 
(c) ^ should scale like ^ — *■ 1 as a — » 0. It turns 
out that criteria (a) is satisfied for mesons, but 
not for baryons; that (b) is marginal; and that 
(c) fails totally (In fact ^ can be large and 
negative !) 

This leads us to the conclusion that the ap- 
plicability of the QCDSR Continuum Model to 
lattice data is questionable since it leads to val- 
ues for So and ^ which are inconsistent with the 
assumptions made. 

It is clear that, apart from taking on non- 
physical values of its fitting parameters, the 
"Cont" Ansatz does reproduce the lattice corre- 
lation function data. Therefore we require similar 
fitting Ansatze as those in the "Cont" case, but 
without the corresponding restrictions on the pa- 
rameters So and ^. 

4. QUARK MODELS 

We present the derivation of the density of 
states, p(s), for hadronic correlators for the non- 
relativistic quark model. Taking the quarks as on- 
shell, there are 3 degrees of freedom available for 
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each quark in the hadron: one per spatial momen- 
tum component. Therefore in the case of mesons 
there are naively 6 degrees of freedom, with 9 
degrees of freedom for baryons. However, there 
are 4 constraint equations corresponding to the 
condition that the hadron's 4-momentum is zero, 
leaving 2 degrees of freedom for the mesons and 5 
degrees of freedom for the baryons. This simple, 
non-relativistic analysis leads to Pnrqm(s) ^ 
for mesons, and Pnrqm[s) ~ for baryons. This 
continuum-like behaviour for p{s) is clearly most 
appropriate for large s where the quarks are ap- 
proximately free. Again we approximate the ac- 
tual density of states by p{s) = 5{s — M) + 
9{s — SQ)pnrqm{s), where the delta-function rep- 
resents the ground state and the non-relativistic 
quark model result is used for the continuum (be- 
ginning at the threshold energy so)- 

Using this definition of p{s) in the spectral 
equation we obtain, for mesons, 

where K is some numerical constant, with a sim- 
ilar expression for baryons. It is very interest- 
ing to note that these expressions are identical 
to those obtained for the Cont model with ruq 
and the condensate set to zero. Since the terms 
0{mq, (: qq ;)) are numerically insignificant, we 
have shown that the non-relativistic quark model 
predicts essentially the same functional form as 
the QCDSR's Continuum Model. Thus we have 
achieved the aim raised in the previous Section of 
a (physically motivated) functional form which 
reproduces the data better than the "2-exp" fits, 
but which doesn't suffer from unphysical parame- 
ter values such as the "Cont" fit. (The relativistic 
quark model is discussed at length in Q.) 

It is natural to wonder if there is any way 
of using the lattice correlation function data to 
pin down the best form of the density of states 
p{s). With this in mind we assume that the 
density of states has the following form: p{s) = 
^S{s — M) + 9{s — so)Ks''^, where n is to be 
determined from the fit. The two-point corre- 
lation function for this Ansatz is defined from 
the standard Laplace transform. Studying the 
values as a function of n we see a clear min- 



imum in the case of the pseudoscalar mesons at 
n = 2. The vector channel's minimum is more 
spread out at around 1 < n < 3. The proton case 
is less clear again; it has a minimum at n ^ 5. 
These results confirm the predictions of the non- 
relativistic quark model. 

5. DISCUSSION & CONCLUSION 

This work began with a study of the QCD 
Sum Rule Continuum Model as applied to lattice 
two-point correlation functions using the method 
first introduced by Leinweber We have ex- 

tended his work by including mesonic states, by 
fitting lattice data at several lattice spacings and 
by using two formulations of the lattice action. 
We have found that the QCD Sum Rule Contin- 
uum Model successfully fits the lattice data, but 
only for unphysical values of its fitting parame- 
ters. This leads us to the conclusion (somewhat 
contrary to that of |^) that this model cannot 
self-consistently fit lattice data. 

We have then searched for a model that repro- 
duces a similar functional form as the QCD Con- 
tinuum Model, but which does not have its re- 
strictions on its parameters. The non-relativistic 
quark model was studied and found to be such a 
theory. We introduced a general fitting function 
which has the non-relativistic quark model as a 
special case, and found evidence that the was 
minimized for the fitting function corresponding 
roughly with the non-relativistic quark model. 

We propose that these fitting functions can be 
used to extract ground state properties in cases 
where the asymptotic state is heavily contami- 
nated by excited states, or where the time sepa- 
ration in the correlators is forced to be small. 
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